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. Abstract 

We analyze the (discrete) spectrum of the semirelativistic "spinless-Salpeter" Hamiltonian 



X 



H = pJm^ + p2 + V{r) , f3>0 



where V{r) is an attractive, spherically symmetric potential in three dimensions. In order to 
5^ ! locate the eigenvalues of H, we extend the "envelope theory," originally formulated only for 
nonrelativistic Schrodinger operators, to the case of Hamiltonians involving the relativistic 
kinetic-energy operator. If V{r) is a convex transformation of the Coulomb potential — 1/r 
and a concave transformation of the harmonic-oscillator potential r^, both upper and lower 
bounds on the discrete eigenvalues of H can be constructed, which may all be expressed in 
the form 

E = min 

r>0 

for suitable values of the numbers P here provided. At the critical point, the relative growth 
to the Coulomb potential h{r) = —1/r must be bounded by dV/dh < 2/3 /n. 
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1 Introduction 

We study the semirelativistic (so-called "spinless-Salpeter" ) Hamiltonian 

H = (3^m? + p2 + V{r) , p>0, (1) 

in which V{r) is a central potential in three dimensions. The eigenvalue equation of this 
operator is called the "spinless Salpeter equation." This equation of motion arises as a well- 
defined standard approximation to the Bethe-Salpeter formalism [Q for the description of 
bound states within a (relativistic) quantum field theory and is arrived at by the following 
simplifying steps: 

1. Eliminate all timelike variables by assuming the Bethe-Salpeter kernel that describes 
the interactions between the bound-state constituents to be static, i.e., instantaneous; 
the result of this reduction step is called the "instantaneous Bethe-Salpeter equation" 
or the "Salpeter equation" 0. 

2. Neglect the spin of the bound-state constituents, assume the Bethe-Salpeter kernel to 
be of convolution type (as is frequently the case), and consider merely positive-energy 
solutions ip, in order to arrive at the so-called "spinless Salpeter equation" Hip = Eip, 
with a Hamiltonian H of the form (||). (For two particles, this form of the Hamiltonian 
H holds only for equal masses m of the bound-state constituents.) 

(For a more detailed account of the reduction of the Bethe-Salpeter equation to the spinless 
Salpeter equation, consult, e.g., the introductory sections of Refs. 0, ^•) This wave equation 
describes the bound states of spin-zero particles (scalar bosons) as well as the spin-averaged 
spectra of the bound states of fermions. 

In this paper we consider potentials which are at the same time convex transformations 
V{r) = g{h{r)) of the Coulomb potential h{r) = — 1/r and concave transformations of the 
harmonic-oscillator potential h{r) = r^. The reason for this is that spectral information is 
known for these two "basis" potentials h{r). Thus the class of potentials is those V{r) that 
have a dual representation 

in which g^^^ is convex ((7(1) > 0) and g^"^^ is concave ((7^^-' < 0). An example of a potential 
in this class is 

V(r) = h C2 In r + csr + c^r'^ , (2) 

r 

where the coefficients {cj} are not negative and are not all zero. Thus tangent lines to the 
transformation function g{h) are of the form ah + b and are either Coulomb potentials lying 
below V, or harmonic-oscillator potentials lying above V. This geometrical idea is the basis 
for our approach to the spectral problem posed by H. We shall consider applications of this 
idea to the (nonrelativistic) Schrodinger problem, the relativistic kinetic-energy operator, 
and the full Salpeter Hamiltonian in Sees. |^, ^, and ^ respectively. We shall show that all 
our upper and lower bounds on the eigenvalues of the semirelativistic Salpeter Hamiltonian 
H of Eq. (IID can be expressed in the compact form 



E ^ min 

r>0 



p2 

/5Wm2 + — + V{r) 
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where P is a constant for each bound, and a sign of approximate equahty is used to indicate 
that, for definite convexity of g{h), the envelope theory yields lower bounds for convex g{h) 
and upper bounds for concave g{h). The main purpose of the present considerations is to 
establish the general envelope formalism in terms of which such bounds can be proved, and 
to determine the appropriate values of P. 

It is fundamental to our method that we first know something about the spectrum of H 
in those cases where V{r) is one of the basis potentials, i.e., the Coulomb and the harmonic 
oscillator. These two spectra are discussed in Sec. |^ below. In Sec. |^ we look at the example 
of the Coulomb-plus-linear potential. 



2 The Coulomb and harmonic-oscillator potentials 
2.1 Scaling behaviour 

Since the two basis potentials are both pure powers, it is helpful first to determine what can 
be learnt about the corresponding eigenvalues by the use of standard scaling arguments. By 
employing a wave function 0(cr) depending on a scale variable c > 0, we find the following 
scaling rule for the eigenvalues corresponding to attractive pure power potentials v sgn(g)r^. 
The Hamiltonian 

H = (3\Jw? + p2 + V sgn(g)r'^ 
has the (energy) eigenvalues E{v,(3,m), where 

E{v,P,m)=f3mE(^-^^,l,iy q > -1 . 

The scaling behaviour described by the above formula allows us to consider the one-particle, 
unit-mass special case m = j3 = 1 initially, that is to say, to work w.l.o.g. with the operator 

H = yl + p^ + f sgn(g)r'' . 



2.2 Coulomb potential 

In the case of the Coulomb potential V{r) = —v/r it is well known that the Hamiltonian 
H has a Friedrichs extension provided the coupling constant v is not too large. Specifically, 
it is necessary in this case that v is smaller than a critical value Vc of the coupling constant: 

2 

V <Vc = - . 

With this restriction, a lower bound to the bottom of the spectrum is provided by Herbst's 
formula 

Eo>^l- {avf , a^^. (3) 

By comparing the spinless Salpeter problem to the corresponding Klein-Gordon equation, 
Martin and Roy have shown that if the coupling constant is further restricted by f < |, 
then an improved lower bound is provided by the expression 



/ 1 + VI - 4t;2 1 
^»^V 2 • "<2- W 
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It turns out that our lower-bound theory has a simpler form when the Coulomb eigenvalue 
bound has the form of Eq. (|^) rather than that of Eq. (^. For this reason, we have derived 
from Eq. (^), by rather elementary methods, a new family of Coulomb bounds. To this end, 
we begin with the ansatz 



>^/l 



[crv] 



and look for conditions under which it becomes true. Since both sides are positive, we may 
square the ansatz and rearrange to yield 



Meanwhile from 
namely. 



Examples are 



we must always satisfy v < ^. This establishes the inequality we'll need. 



a 



a 



lav 



1 

"^2 



V < 



1 1 
-■=2 



(5) 



1.145898 



.<3; 



^2 = 8-4^3 ^ 1.071797 , v<\. 

~ 4 

All these (lower) bounds are slightly weaker than the Martin-Roy bound but above the 
Herbst bound (|^). We note that these functions of the coupling constant v are all monotone 
and concave. 



2.3 Harmonic-oscillator potential 

In the case of the harmonic-oscillator potential, i.e., V{r) = vr"^, much more is known |^ 
In momentum-space representation the operator p becomes a c-variable and thus, from the 
spectral point of view, the Hamiltonian 

H = y^l + p2 + vr^ 

is equivalent to the Schrodinger operator 

H = -vA + Vl +r2 . (6) 

Since the potential in this operator increases without bound, we know ||^ that the spectrum 
of this operator is entirely discrete. We call its eigenvalues £ne{v), n = 1, 2, ...,£ = 0, 1, ... , 
where n counts the radial states in each angular-momentum subspace labelled by i. In what 
follows we shall either approximate the eigenvalues Sniiv) analytically or presume that they 
are known numerically. The concavity of nonrelativistic Schrodinger energy eigenvalues has 
been discussed in Refs. |jlO|, |l| . Theorem 2 of Ref . [|12] establishes concavity for the ground 
state; the same proof can be applied to states which are (1) in the subspace corresponding to 
angular momentum i and (2) orthogonal to the first n — 1 exact energy eigenstates in this 
subspace. This establishes concavity also for all the higher Schrodinger energy eigenvalues. 
Thus the eigenvalues £ne{v), regarded as functions of the coupling parameter v, are concave. 
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2.4 The spectral comparison theorem 

For the class of interaction potentials given by (Q) with the coefficient of the Coulombic term 
satisfying the constraint 

limr^V'ir) < — , 

the semirelativistic Salpeter Hamiltonian H is bounded below and is essentially self-adjoint 
1^. Consequently, the discrete spectrum of H is characterized variationally |^ and it follows 
immediately from this that, if we compare two such Hamiltonians H having the potentials 
V'^^^r) and V^'^\r), respectively, and we know that V^^\r) < V(2)(r),then we may conclude 
that the corresponding discrete eigenvalues Ene satisfy the inequalities E^^^ < E^) . We shall 
refer to this fundamental result as the "spectral comparison theorem." In the more common 
case of nonrelativistic dynamics, i.e., for a (nonrelativistic) kinetic term of the form /?p^/2m 
in the Hamiltonian H, a constraint similar to the above would hold for the coefficient of a 
possible additional (attractive) — term in the potential V{r). 



3 Envelope representations for Schrodinger operators 

We distinguish a potential V{r) = vf{r) from its shape f{r), where the positive parameter v 



is often called the "coupling constant." The idea behind envelope representations |jT2|, |T3| is 
suggested by the question: if one potential /(r) can be written as a smooth transformation 
/(r) = g{h{r)) of another potential h{r), what spectral relationship might this induce? We 
consider potential shapes that support at least one discrete eigenvalue for sufficiently large 
values of the coupling v and suppose for the sake of definiteness that the lowest eigenvalue of 
—A-\-vh{r) is given by H{v) and that of —A+vf{r) by F{v). If the transformation function 
g{h) is smooth, then each tangent to g is an affine transformation of the "envelope basis" h 
of the form f^^\r) = a{t)h{r)+b(t), where r = t is the point of contact. The coefficients a{t) 
and b{t) are obtained by demanding that the "tangential potential" /^*^ (r) and its derivative 
agree with /(r) at the point of contact r = t. Thus we have 



a{t) = {|| , b{t) = fit) - a{t)h{t) . 

The corresponding geometrical configuration is illustrated in Fig. |1| in which the potential / 
is chosen to be the Coulomb-plus-linear potential, /(r) = — 1/r+r, and the envelope basis h 
is, for the upper family, the harmonic-oscillator potential h{r) = and, for the lower family, 
the Coulomb potential h{r) = — 1/r. The spectral function F^^^ for the tangential potential 
yW(r) = a{t)h{r) + h{t) is given by F^^\v) = H{va{t)) + vh{t). If the transformation g{h) 
has definite convexity, say g"{h) > 0, then each tangential potential f^^\r) lies beneath /(r) 
and, as a consequence of the spectral comparison theorem, we know that each corresponding 
tangential spectral function F^^\v), and the envelope of this set, lie beneath F{v). Similarly, 
in the case where g is concave, i.e., g"{h) < 0, we obtain upper bounds to F{v). These purely 
geometrical arguments, depending on the spectral comparison theorem, extend easily to the 
excited states of the problem under consideration. The spectral curves corresponding to the 
envelope representations for the potential in Fig. |I] are shown in Fig. ^ for the excited state 
(n, £) = (2, 4). For comparison the exact curve E = F{v) is also shown in Fig. this curve 
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will be close to the Coulomb envelope for large v and to the oscillator envelope for small v. 
Of course, the envelopes of which we speak still have to be determined explicitly. Extensions 
of this idea to completely new problems, such as simultaneous transformations of each of a 
number of potential terms [Q, the Dirac equation [^, or the spinless-Salpeter problem 
of the present paper, are best formulated initially with the basic argument outlined above. 



f(r) 




12 3 



r 

Figure 1: The "Coulomb-plus-linear" potential shape /(r) = —a/r + br represented as the 
envelope curve of two distinct families of tangential potentials of the form ah{r) +j3. In the 
upper family h{r) = is the harmonic-oscillator potential; in the lower family h{r) = — 1/r 
is the Coulomb potential. The adopted values of the relevant physical parameters a and b are 
a = 0.2 and b = 0.5. 



In the 1-term case the question arises as to whether there is a simple way of determining 
the envelopes of the families of upper and lower spectral functions. One effective solution of 
this problem is by the use of "kinetic potentials" which were introduced |T2|, |13| precisely for 
this purpose. The idea is as follows. To each spectral function F^iiv) there is a corresponding 
"kinetic potential" (that is, a minimum mean is o-kinetic potential) fni{,s). The relationship 
between F„£ and fne is invertible and is essentially that of a Legendre transformation |17]: we 
can prove in general that F is concave, f{s) is convex, and 

F"{v)f"is) = -1 . 
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V 

Figure 2: The spectral approximation corresponding to Fig. |l|. Each "tangential" potential 
y(t)(^r) = ah{r)+P generates a corresponding tangential energy curve 

The envelopes of these spectral families generate upper and lower bounds to the exact curve 
E = F{v), shown here for the case {n,i) = (2,4). 



The explicit transformation formulas are as follows: 

fUs) = F:,,iv), s = FUv)-vFM (7) 

and 

= fni{s) - Sf'^^{s) , - = -fU{s) . (8) 

An a-priori definition of the ground-state kinetic potential /io(s) = f{s) is given by 

f{s) = inf f^) , 

(V.,t/.)=i 

(i/',-Ai/')=s 

where T> C L^(IR^) is the domain of the Hamiltonian. The definition for the excited states is 
a little more complicated [0 and, in view of (0) and (H), will not be needed in what follows. 

What is crucial is that the spectral functions, either exact or approximate, are recovered 
from the corresponding kinetic potentials by a minimization over the kinetic-energy variable 
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s. In this way the total minimization required by the minimum-maximum principle |TI 
has been divided into two steps: the first is constrained by [ip, — A^/)) = s and the second is a 
minimization over s. We have in all cases: 



FnAv) 



mm 

s>0 



S + Vfnl{s) 



Another form of this expression is possible for the kinetic potential is monotone and allows 
us to change variables (s — > r) by /(r) = fne{s). Thus we have 



Fniiv) 



mm 

r>0 



if) 



(9) 



The two corresponding expressions of the envelope approximation then become 



fni{s) ^ g{Ki{s)) , or ir, 



{h) 

n£ 



The second form (|]) of the explicit expression for F„^(t>) isolates the potential shape / itself 
and leads to an inversion sequence |]TB[ which reconstructs the potential from a single given 



spectral function; but this is another story |19, 20 



It is useful here to provide the formulas for the kinetic potentials / corresponding to pure 
power-law potentials V{r) = v sgn{q)r'^. Elementary scaling arguments for the Hamiltonian 



show that the eigenvalues satisfy 



H = — A + V sgn{q)r'^ 



Fniiv) = v—^FMl) 



where Fne{l) = Ene{q) are the eigenvalues of H with coupling v = 1, i.e., of — A + sgn(g)r'^. 
From Eq. (0) we immediately find that the kinetic potentials fne{s) for these potentials V{r) 
are given by 



fn 



qEnM) 



2 + q 



Meanwhile the corresponding i^-functions all have the same simple form |2T 

2 



2 + q 
. ha? 

'PnM) 



S 2 



(10) 



where the P numbers are given by 



Pni{q) = \Eni{q)\ 



2+q 
2q 



2 + 



2 + q 



qy^O 



Consequently, the power-law kinetic potentials (|T0D may be expressed in the simple form 



fn 



sgn((?) 



q^O. 



Some of the eigenvalues E, and thus the corresponding P numbers, are known exactly from 
elementary quantum mechanics. From the known eigenvalues E for the Coulomb potential. 
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En£{—^) = ~[2('^ + ^)] ^, and the harmonic-oscillator potential, Eni{2) 
immediately obtain the corresponding P numbers: 



4n + 2£- 1, 



we 



n + 



Pn^(2) 



2n + ^- 



The case g = corresponds exactly to the logarithmic potential ; the P numbers for 

the lowest states of the logarithmic and the linear potentials may be found in Table ||. 



Table 1: Numerical values for the P numbers for the logarithmic potential (g = 0) and the 
linear potential (g = 1) used in the Schrodinger eigenvalue formula ([TT|). 



Ti 








1 





1.21867 


1.37608 


2 





2.72065 


3.18131 


3 





4.23356 


4.99255 


4 





5.74962 


6.80514 


5 





7.26708 


8.61823 


1 


1 


2.21348 


2.37192 


2 


1 


3.68538 


4.15501 


3 


1 


5.17774 


5.95300 


4 


1 


6.67936 


7.75701 


5 


1 


8.18607 


9.56408 


1 


2 


3.21149 


3.37018 


2 


2 


4.66860 


5.14135 


3 


2 


6.14672 


6.92911 


4 


2 


7.63639 


8.72515 


5 


2 


9.13319 


10.52596 


1 


3 


4.21044 


4.36923 


2 


3 


5.65879 


6.13298 


3 


3 


7.12686 


7.91304 


4 


3 


8.60714 


9.70236 


5 


3 


10.09555 


11.49748 


1 


4 


5.20980 


5.36863 


2 


4 


6.65235 


7.12732 


3 


4 


8.11305 


8.90148 


4 


4 


9.58587 


10.68521 


5 


4 


11.06725 


12.47532 


1 


5 


6.20936 


6.36822 


2 


5 


7.64780 


8.12324 


3 


5 


9.10288 


9.89276 


4 


5 


10.56970 


11.67183 


5 


5 


12.04517 


13.45756 
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In summary, if the potential V{r) is a smooth transformation V{r) 
pure power-law potential sgn(g)r'^, then the eigenvalues of 

H = -A + V{r) 

are given approximately by the expression 



g{sgn{q)r'^) of the 



Eni ^ min 



r>0 



+ Vir) 



Here, a sign of approximate equality is used to indicate that, for a definite convexity of g{h), 

Eq. (|1T]) yields lower bounds for convex g {g" > 0) and upper bounds for concave g {g" < 0). 

The numbers Pni{<l) can be derived from the eigenvalues of the operator —A + sgn(g)r''. 
The lower bounds derived in the framework of envelope theory can be improved by use of 

the refined comparison theorems of Ref. , which allow comparison potentials to intersect; 

a detailed study of the latter approach is, however, beyond the scope of the present analysis. 
As an immediate application we consider the (nonrelativistic) Schrodinger Hamiltonian 
for the semirelativistic spinless-Salpeter harmonic-oscillator problem (|l]). Here we have 



H = -vA + V{r) 



with 



V{r) = (3^m? + r2 ; 

this potential is a convex transformation of a linear potential and a concave transformation 



of a harmonic-oscillator potential. We conclude therefore from Eq. (|TT]) (see also Ref. [p5[1): 



mm 

r>0 



< min 

T->0 



+ j3\/m? + 



(12) 



where the numbers PniiX) given in Table |T] and Pne{2) = 2n + i—^. By a simple change 
of variables, r r' = P/r, we are able to recast the inequalities (p!^) into the "preferred" 
form 



n£[V, 



^ mm 

r>0 



/3\m^ + ^ + Vir) 



V{r) 



vr 



in which the function to be minimized is simply the spinless-Salpeter Hamiltonian with the 
momentum operator |p| replaced by P/r; the P numbers yielding upper and lower bounds 
are as in Eq. (0). Interestingly, the upper and lower bounds in Eq. (|^) are equivalent to the 
corresponding bounds obtained in Ref. p|; however, these earlier specific bounds were not 
derived as part of the general envelope theory. If we approximate the square root from above 
again, by using the elementary inequality 



vm^ + r^ < m H , 

2m 



we obtain from (|12|) the weaker upper bound 



Sneiv) < /?m + y'^(4n + 2£ - 1) , 
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which is identical to that given by the general "Schrodinger upper bound" |, ^ obtained 
by initially approximating the relativistic kinetic-energy operator above by 

The latter upper bound, and an improvement on it, will be discussed in the next section. 



Envelope approximations for the relativistic kinetic 
energy 



The relativistic kinetic-energy operator T = (5\/m^ + is a concave transformation of p^. 
Thus "tangent lines" to this operator all are of the form ap^ -|- h and each one generates a 
Schrodinger operator that provides an upper bound to T. In a given application with given 
parameter values, one can search for the best such upper bound. By elementary analysis we 
can establish the operator inequality 



A 



E = (3Jm? + p2 + V{r) + ^+ _ + V{r) , 



(13) 



where /i = \ w? and |p| = pi is the "point of contact" of the tangent line with the 



square-root function. The inequality (p^ ) is identical to that obtained by employing the 
inequality (T — //)^ > 0. Optimization over yU for the Coulomb case V{r) = —v/r recovers 
the explicit upper-bound formula of Ref. 



Eniiv) < mP^ 

\ \p[n + t 

In the case of the harmonic-oscillator potential, V{r) = f r^, we obtain the upper bounds 



Ene{v) < min 

/i>0 



|;(4n + 2«-l) + fL+- 



(Brief reviews of analytical upper bounds on the energy eigenvalues of the spinless Salpeter 
equation derived by combining operator inequalities with the minimum-maximum principle 
may be found in Refs. 0, ^, |27[|.) 



The strategy of the present section is to regard the relativistic kinetic-energy operator T 
as a concave function of p^, so that "tangent lines" generate "upper" Schrodinger operators. 
This general approach leads to the same upper bounds as those of Martin pB| who used the 
particular square-root form of the relativistic kinetic energy to construct an operator whose 
positivity yields the bounds. 



5 Envelope approximations for Salpeter Hamiltonians 
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5.1 The principal envelope formula 

Let us now turn to our main topic and consider the spinless-Salpeter Hamiltonian of Eq. 

H = /3^m2 + p2 + V{r) , 

and its eigenvalues E. We shall assume that the potential V{r) is a smooth transformation 
V{r) = g{h{r)) of another potential h{r) and that g has definite convexity so that we obtain 
bounds to the energy eigenvalues E. We suppose that the "basis" potential h{r) generates a 
"tangential" Salpeter problem 

H = f3\Jm? + p2 + vh{r) , 

for which the eigenvalues e(f ), or bounds to them, are known. We shall follow here as closely 
as possible the development in Sec. |^ for the corresponding Schrodinger problem. 

First of all, we recall that the approximations or bounds to the energy eigenvalues of the 
relativistic Coulomb and harmonic-oscillator problems we shall eventually use from Sees. |^ 
and ^, when regarded as functions of the coupling f , are all concave. Furthermore, it is easy 
to convince oneself that all the (unknown) energy functions e{v) of the "tangential" Salpeter 
problem are concave, that is, e"{v) < 0. Suppose that the exact eigenvalue and (normalized) 
eigenvector for the problem posed by the "tangential" Hamiltonian 



Ti = l3\lm? + p2 + vh{r) 

are e{y) and ip{y, r). Then, by differentiating the expectation value (■?/', TLiIj) with respect to 
the coupling t>, we find 

e'{v) = {tpy hip) . 

If we now apply ip{v,r) as a trial vector to estimate the energy of the operator 



I3\jm? + p2 + uh{r-) , 

in which v has been replaced by u, we obtain an upper bound to the corresponding energy 
function e{u) which may be written in the form 

e{u) < e{v) + {u — v)e'{v) . 

This inequality tells us that the function e{u) lies beneath its tangents; that is to say, e{u) is 
concave. Convexity properties of the energy functions of the corresponding (nonrelativistic) 



Schrodinger problem have been investigated in Refs. |1T0| , |12| pT| . 

Next, in order to prove the main result of this section, the "principal envelope formula," 
we begin by using an envelope representation for the potential V{r) in the Hamiltonian (|I|) 
and then demonstrate that all the spectral formulas that follow possess a certain structure. 
Finally, as an application, we specialize to the case of pure power-law "basis" potentials h{r) 
and, more particularly, to the Coulomb potential and the harmonic-oscillator potential for 
which, at this time, we have spectral information [cf. the discussions in Sees. |2.2| and 



and the exact bounds (|T2D on the energy levels Sniiv) of the relativistic harmonic oscillator]. 
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The tangential potentials we shall employ have the form V'^^\r) = a{t)h{r)+b{t), where, 
as in the Schrodinger case, the coefficients a{t) and b{t) are given by 

a{t) = ^ = 9'iHt)) , bit) = Vit) - a{t)h{t) = g{h{t)) - g'{h{t))h{t) , 

and r = t is the point of contact of the potential V{r) and its tangent V'^^\r). If, for the sake 
of definiteness, we assume that V = g{h) with g concave (i.e., g" < 0), we obtain a family of 
upper bounds given by 

E < 6{t) = e(a(t)) + b{t) . 
The best of these is given by optimizing over t: 

E < 6{i) = e(a(t)) + b{i) , 

where t, the value of t which optimizes these bounds, is to be determined as the solution of 

e'ig'ihii))) = h{i) . 



In the spirit of the Legendre transformation we now consider another problem which has 
the same solution; this second problem is the one that provides us with our basic eigenvalue 
formula. We consider 

S = nnn[e(f ) — ve'{v) + g{e'{v))] , 
which is well defined since e{v) is concave. The solution has the critical point 

V = g'{e'{v)) . 

If we now apply the correspondence h{i) = e'{v), it follows that the critical point v becomes 

V = g\h{i)) , 

and the tangential-potential coefficients a and b become 

a{i)=g\e\v)) = v , b{i) = g{e'{v)) - ve'{v) , v = v. (14) 
Meanwhile the original critical (energy) value is given by 

e{i) = e{a{i)) + b{t) = e{v) - ve'{v) + g{e'{v)) , v = v . 

Thus we conclude that the spectral approximation obtained by envelope methods is given 
by the following "principal envelope formula:" 

E^S = min[e(f ) - ve'{v) + g{e'{v))] . (15) 

t)>0 

If g is concave (that is, g" < 0), then E < S; if g is convex (that is, g" > 0), then E > S. 
From the above considerations it follows immediately that, if the exact energy function e{v) 
corresponding to the basis potential h is not available, then, for g{h) concave, concave upper 
approximations eu(f ) > e{v) or, for g{h) convex, concave lower approximations ei{v) < e{v) 
may be used instead of the exact energy function e{v) in the principal envelope formula (|T3p. 
Then all the lower tangents will lie even lower and all the upper tangents will lie even higher. 
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If g is convex, we obtain a lower bound; if g is concave, we obtain an upper bound; because of 
the concavity of e{v), this extremum is a minimum in both cases. If we wish to use numerical 
solutions to the "basis" problem (generated by h{r)), or if a completely new energy-bound 
expression becomes available, the principal envelope formula (^) is the relation that would 
at first be used. 

Interestingly, in the formula (0) the tangential-potential apparatus is no longer evident; 
only the correct convexity is required. As in the Schrodinger case []T^ , once we have the basic 
result, the reformulation in terms of "kinetic potentials" is often useful: the kinetic potential 
h{s) corresponding to some basis potential h{r) is given by the Legendre transformation fl^ 

h{s) = e'{v) , s = e{v) — ve'{v) . 

Meanwhile the envelope approximation has the kinetic-potential expression V{s) ~ g{h{s)). 

For both the Coulomb lower bounds (^) or (^ and the harmonic-oscillator upper bounds 
(|1^) which we have at present, we may express our general results in a special common form 
which will now be derived. 



5.2 The Coulomb lower bound 



We consider first the Coulomb lower bound in which we assume that the potential V{r) is a 
convex transformation V{r) = g[h{r)) of the Coulomb potential h{r) = —1/r. According to 
Sec. in this case all the "lower" ei{v) have been arranged — with the parameters /3 and m 
returned — in the form 



ei(f ) = Pm 



\ 



From this it follows by elementary algebra that if we define a new optimization variable r by 

e[{v) = h{r) = —1/r, we have 



ei(f) — ve[{v) = (3\ m'^ + 



p2 



P 



1 

cr 



Consequently, the lower bound on the energy eigenvalues E of the spinless Salpeter equation 
becomes 



E > min 

r>0 



V < (3vp . 



(16) 



Here the boundary value of the Coulomb coupling v is given, when simply determined by 
the requirement of boundedness from below of the operator (|I]), by the critical coupling Vc, 



Vp = Vc 



2 

TT 



and, when arising from the region of validity of our Coulomb-like family of lower bounds (^, 
via P = l/o", by 

vp = PVT^<^. (17) 

Some {P, Vp} pairs may be found in Table ^; others can be easily generated from the upper 
bound on the coupling v given in Eq. {^7^. The meaning of the Coulomb- coupling constraint 
is a{t) < f3vp, where a is the coefficient in the tangential Coulomb potential given by (p!4D. 
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Table 2: Explicit values of some {P, vp} pairs, obtained via the equality P = l/cx either from 
the Herbst lower bound or the expression (|^) for our new lower bounds on the spectrum 
of the spinless relativistic Coulomb problem (in three spatial dimensions) . 



P vp 

2 2 

vr vr 

1 1 

V 9 - 3^5 3 

1 1 

2^2-^3 ^ 



As a rather trivial consistency check of our formalism, the Coulomb lower energy bound 
of Eq. (|16D may be applied to the Coulomb potential V{r) = —v/r in order to re-derive, for 
P = 2/7T, the Herbst formula (P) — which is nothing else but the starting point of the present 
"lower-bound" considerations. 

5.3 The harmonic-oscillator upper bounds 

Next, let us turn to the harmonic-oscillator upper bounds. Our main assumption is here that 
V{r) = g{r'^)i with g" < 0. In this case the only difficulty is that the basis problem h{r) = 
is equivalent to a Schrodinger problem whose solution Sne{v) is not known exactly. Following 



the discussion of suitable bounds after the proof of the principal envelope formula, Eq. (|T5[), 
let us call the upper bound provided by Eq. (p!^ ) eu(f), and let us introduce the shorthand 
notation Pne{2) = 2n+£— | = P. Then we have the following parametric equations for eu(f ): 

en{v) =v— + pVm'^ + , v = ——^=== , e'^{v) = — . 

2P^\/m? + r'^ 

By substituting these expressions into the fundamental envelope formula (p!5D we obtain the 
following upper bound on all the eigenvalues of the spinless-Salpeter problem with potential 
V{r) = g{r'^) and g" < 0: 



Ene < min 

r>0 



P = Pne{2) = 2n + i-^. (18) 



6 The Coulomb-plus-linear (or "funnel") potential 

In order to illustrate the above results by a physically motivated example, let us apply these 
considerations to the Coulomb-plus-linear or (in view of its shape) "funnel" potential 

V{r) = + Car , Ci > , Cs > . 
r 
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(This potential provides a reasonable overall description of the strong interactions of quarks 
in hadrons. For the phenomenological description of hadrons in terms of both nonrelativistic 
and semirelativistic potential models, see, e.g., Refs. ||2^, 0.) By choosing as basis potential 
the Coulomb potential h{r) = — 1/r, we may write V{r) = g{h{r)) with 

g{h) = ci/i - ^ , 

which is clearly a convex function of /i < 0: g" > 0. Thus the convexity condition is satisfied. 
However, we are not free to choose the coupling constants Ci and C2 as large as we please. It is 
immediately obvious that, for a particular {P, vp} pair, we must in any case have ci < jSvp. 
For the full problem the coefficient C2 of the linear term will also be involved. The coupling v 
we are concerned about is given by (p!^. We have 

V = g[e [v]) = — , _ -- = ci + — = ci + Car . 




From this we obtain, for given values of the parameters m and (3 and for a given {P, vp} pair, 
as a sufficient condition for v < (3vp the "Coulomb coupling constant constraint" on the two 
coupling strengths Ci and C2 in the funnel potential: 

p2 /p2 \ 

c, + -(^-^-ljc2</?.p. (19) 

In the case {P = 1/V2, vp = 1/2} and P = m = 1 this condition reduces to Ci + |c2 < |. For 
Herbst's lower bound (|^), i.e., P = vp = Vc = 2/7r, this constraint clearly yields ci < Pvp. 
There is no escaping this feature of all energy bounds involving the Coulomb potential: the 
constraint derives from the fundamental observation that the Coulomb coupling v must not 
be too large, so that the (relativistic) kinetic energy is able to counterbalance the Coulomb 
potential in order to maintain the Hamiltonian (|I]) with V{r) = —v/r bounded from below. 

For example, if we seek the largest allowed value of the parameter P by solving Eqs. ([T7|) 
and ( [I9| ) together, we find that this largest P is given by 



C2 sin"^ t 2 



cos^ t{l3 sin t cos t — ci) 



m% P = sint. (20) 



For the Coulomb-plus-linear potential V{r) = — ci/r -|- C2r under consideration. Fig. |^ 
shows the lower and upper bounds on the lowest energy eigenvalue E of the spinless Salpeter 
equation, given by the envelopes of the lower and upper families of tangential energy curves 
(|1^) and (p!8D, as functions E{m) of the mass m entering in the semirelativistic Hamiltonian. 



In the case of the Coulomb lower bound (16), we have employed for each m the best possible 



P{m) provided by (pOl). As m 0, the "basis" Coulomb problem H = [3\Jm? + -p^—v/r has 
energy e(m) — > 0; thus the Coulomb lower bound for a non-Coulomb problem becomes very 
weak for small values of m. Of course, Eq. (p!8| ) provides us with rigorous upper bounds for 
every energy level. 

In order to get an idea of the location of the exact energy eigenvalues E, Fig. ^ also shows 
the ground-state energy curve E{m) obtained by the Rayleigh-Ritz variational technique 
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Figure 3: Lower bounds (L), according to (|T6D, and upper bounds (U), according to (|T^), on 
the energy eigenvalue E of the ground state [{n, i) = (1, 0)] of the spinless Salpeter equation 
with a Coulomb-plus-hnear potential V{r) = —ci/r+C2r, for /5 = 1, Ci = 0.1, and C2 = 0.25. 



The lower bound is given by the general result ([TBD with the "best" P{m) provided by (pOf). 
In order to satisfy the Coulomb coupling constraint (PDj), the mass m must fulfil m > Vb/A. 
For comparison, a (very accurate) Rayleigh-Ritz variational upper bound E is depicted too. 



with the Laguerre basis states for the trial space defined in Ref. [0. Strictly speaking, this 
energy curve represents only an upper bound to the precise eigenvalue E. However, from the 
findings of Ref. the deviations of these Laguerre bounds from the exact eigenvalues may 
be estimated, for the superposition of 25 basis functions used here, to be of the order of 1 %. 

Figure ^ shows, for a Coulomb-plus-linear potential of the form V{r) = v{—a/r+br), the 
lower and upper bounds on the lowest energy eigenvalue E of the spinless Salpeter equation, 
given by the envelopes of the lower and upper families of tangential energy curves (|I6]) and 
(|l^), as functions E{v) of the "overall" coupling parameter v which multiplies the potential 
shape —a/r + br. Again we compare these bounds with the ground-state energy curve E{v), 
obtained by the Rayleigh-Ritz variational technique with the Laguerre basis states |]3T| . 
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Figure 4: Lower bounds (L), according to (p!6D, and upper bounds (U), according to (0), on 
the energy eigenvalue E of the ground state [{n, i) = (1, 0)] of the spinless Salpeter equation 
with Coulomb-plus-hnear potential V{r) = v{—a/r + br), for a = 0.2, b = 0.5, m = (3 = 1. 
The lower bound is given by the general result (|16D with P(m = 1) computed from Eq. (pOD, 
evaluated for mass m = 1 and the funnel-potential coupling strengths Ci = va and C2 = vh. 
For comparison, a (very accurate) Rayleigh-Ritz variational upper bound E is depicted too. 



7 Summary and conclusion 

In this analysis we have studied the discrete spectrum of semirelativistic "spinless-Salpeter" 
Hamiltonians H, defined in Eq. (|1|) , by an approach which is based principally on convexity. 
We have at our disposal very definite information concerning, on the one hand, the bottom 
of the spectrum of H for the Coulomb potential, h{r) = — 1/r, and, on the other hand, the 
entire spectrum of H for the harmonic-oscillator potential, h{r) = . The class of potentials 
that are at the same time a convex transformation of —1/r and a concave transformation of 
includes, for example, arbitrary linear combinations of Coulomb, logarithmic, linear, and 
harmonic-oscillator terms. In order to obtain information about the eigenvalues E of H for 
arbitrary members within this class of potentials, we have extended the — for nonrelativistic 
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Schrodinger operators well-established — formalism of envelope theory to Hamiltonians with 
relativistic kinetic energies. The envelope technique applied here takes advantage of the fact 
that all "tangent lines" to the interaction potential V{r) = g{h{r)) in H are potentials of the 
form ah{r) + b, and that, by convexity and the comparison theorem recalled in Subsec. |2.4| , 
the energy eigenvalues corresponding to these "tangent" potentials provide rigorous bounds 
to the unknown exact eigenvalues E of H. If e{v) denotes the energy function — or a suitable 
bound to it — corresponding to the problem posed by a "basis" potential vh{r), where f is a 
positive coupling parameter, the envelopes of upper and lower families of energy curves may 
be found with the help of the "principal envelope formula" 



Here, a sign of approximate equality is used to indicate that, for a definite convexity of g{h), 
the envelope theory yields lower bounds for convex g{h) and upper bounds for concave g{h). 
With the above principal envelope formula at hand, all new spectral pairs {h{r), e{v)} which 
may become available at some future time can immediately be used to enrich our collection 
of energy bounds. If the basis potential h{r) is a pure power, these bounds can be written as 



where the numbers Pni are obtained from the corresponding underlying basis problems. The 
power of this technique is illustrated, in Sec. ^ by our application to the "funnel" potential, 
V{r) = —ci/r + C2r. For this problem, we have employed both the semirelativistic Coulomb 
and harmonic-oscillator problems to calculate, respectively, lower and upper bounds on the 
energy eigenvalues of the spinless Salpeter equation. 

We expect that such results would provide bounds on the energy eigenvalues for general 
theoretical discussions, or be used as guides for more tightly focussed analytic or numerical 
studies of the spectra of semirelativistic "spinless-Salpeter" Hamiltonians. 
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